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Abstract 

In this article we investigate the linear-chain spin-half J\—Ji model by using high-order coupled 
cluster method (CCM) calculations. We employ three model states, namely, a nearest-neighbour 
(n.n.) Neel model state in which neighbouring spins are antiparallel, a next-nearest-neighbour 
(n.n.n.) Neel model state in which next-neighbouring spins are antiparallel, and finally a type of 
"double spiral" model state with two pitch angles. For J2/J1 < 5, we find that the n.n. Neel model 
state produces the lowest energies. For J2/J1 > \, we find that the stable states for the quantum 
system are those for the "traditional" spiral state in which the two pitch angles are identical and 
the collinear n.n.n. Neel model state. As seen previously, we show that we are able to reproduce 
exactly the dimerised ground (ket) state at the Majumdar- Ghosh point (J2/J1 = \) using the n.n. 
Neel model state. We show that the onset of the dimerised phase is indicated by a bifurcation 
of the nearest-neighbour ket- and bra-state correlation coefficients for the nearest-neighbour Neel 
model state. Furthermore, we show that the n.n.n. Neel model state can also reproduce exactly the 
dimerized ground (ket) state at the Majumdar-Ghosh point (J2/J1 = \)- The ground-state energies 
for the n.n.n. Neel model state are shown to be much lower than those of the "spiral" model state 
for a finite region above J2/J1 = \- Indeed, we show that ground-state energies for the collinear 
n.n. and n.n.n. model states are in good agreement with the results of exact diagonalisations of 
finite-length chains across this entire regime for J\ > 0. We produce results for the dimer order 
parameter for J2/J1 > \ by using the n.n.n. Neel model state and these results are in reasonable 
agreement with previous results of DMRG. Finally, results for the spin-spin correlation functions 
for the n.n.n. Neel model state are presented. It is shown that the correlation length increases 
with increasing J2/J1 above J2/J1 = \ using this model state and that shift in the peak of the 
structure function occurs from q , = 7rtog = 7r/2as one increases J2/ J\- Although these results are 
not conclusive, they are intriguing because they suggest that a ground state exhibiting collinear 
order might occur for some region for J2/J1 > \- 
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(a) Spiral 



(b) n.n. Neel 



(c) n.n.n. Neel 




FIG. 1: Model states for the one-dimensional spin-half J\—Ji model, (a) Spiral model states with 
two angles; (b) collinear nearest-neighbour (n.n.) Neel model state; and, (c) collinear next-nearest- 
neighbour (n.n.n.) Neel model state 



I. INTRODUCTION 



The formation of dimer- and plaquette-ordered singlet ground states (so-called valence- 
bond crystal (VBC) states) is an interesting and important topic in quantum spin systems. 
Often, the formation of enhanced dimer or plaquette correlations is driven by frustration, 
which can increase quantum fluctuations and which may result in such gapped rotationally- 
invariant quantum paramagnetic states. An example for the appearance of such exact VBC 
product eigenstates are the spin-half J1-J2 model on the linear chain [1-12] at the point 
J2/ J\ = \ (the so-called Majumdar-Ghosh point), which is an example of spontaneous 
symmetry breaking. 

The Hamiltonian for this spin-half model has nearest-neighbour (n.n.) bonds of strength 
Ji and next-nearest-neighbour (n.n.n.) bonds of strength J2 is given by 



H 



■h \ - 



i+pi ~r Si ' S i+P2 ) 

i,p2 



(1) 



where the index i runs over sites on the lattice, p\ runs over all nearest-neighbours to i, 
and P2 runs over all next-nearest-neighbours to %. Henceforth we put | J x | = 1 and consider 
■h > 0. 



The ground-state properties of this system have been studied using methods such as 
exact diagonalizations [2,7], DMRG [3-5,9], CCM [8-11], and field-theoretical approaches [5] 
(see Refs. [5,6] for a general review). At J2/J1 = we have the unfrustrated Heisenberg 
antiferro magnet, where the exact solution is provided by the Bethe Ansatz. The ground state 
is gapless and the spin-spin correlation function (sj • Sj) decays slowly to zero according to 
a power-law, i.e. no true Neel-like long-range order is observed. In the region J2/J1 > 
the nearest-neighbour (Ji) and next-nearest-neighbour interactions (J2) compete, thus 
leading to frustration. At J2/J1 = 0.2411(1) the model exhibits a transition to a two-fold 
degenerate gapped dimerized phase with an exponential decay of the correlation function 
(Sj • Sj) [2-6]. This state breaks the translational lattice symmetry. At the Majumdar- 
Ghosh point J2/J1 = \ there are two degenerate simple exact dimer-singlet product ground 
states corresponding to the dimerized product state for the Hamiltonian of Eq. ([I]) [1]. The 
ground state for J2/J1 > \ is not completely well- understood, although it is known that 
gap [5] persists for some range above J2/J1 = | and that incommensurate spiral correlation 
occur. Furthermore, DMRG studies suggest that a Liftschitz point occurs at J2/J1 = 0.52 
and that the dimer order parameter peaks at J2/J1 = 0.58. In the limit J2/J1 — > 00, the 
system decouples into two Bethe chains and so the system is gapless. 

II. METHOD 

A. The CCM Formalism 

A description of the underlying methodology for the CCM is given in Refs. [13-21]. We 
note here that the ket and bra ground-state energy eigenvectors, and of a general 
many-body system described by a Hamiltonian H 

H\$) = E g \V); (<I>\H = E g {*\ , (2) 
are parametrised within the single-reference CCM as follows: 

|vf} = e 5 |$> ; S = J2SiCt, 
(§| = ($|^e- 5 ; S = l + J2SiC7. (3) 

Note that |$) is the normalised single model or reference state. The correlation operator S 
is thus a linked-cluster operator that is decomposed in Eq. (j3J) in terms of a complete set 
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of mutually commuting multi-spin and multi-configurational creations operators Cj with 
respect to the model state. The hermitian adjoints, Cj = (Cfy are the corresponding 
destruction operators, ($|C/~ = = C*7|$), and we explicitly define Cq = 1, the identity 
operator. The label / is thus a set-index comprising a set of single-particle labels in some 
suitable single-particle basis defined via |$). We note that the normalisation conditions 
( < &| V I / ) = ($|$) = 1 = (^l^) follow from Eq. ([3]). The ground ket- and bra-state equations 
are given in terms of H = (^f\H\^f) as 



In order to solve for the one-dimensional J1-J2 model considered here, we make approx- 
imations in both 5" and 5". The three most commonly employed approximation schemes 
previously utilised have been: (1) the SUBn scheme, in which all correlations involving only 
n or fewer spins are retained, but no further restriction is made concerning their spatial 
separation on the lattice; (2) the SUBn-m sub-approximation, in which all SUBn correla- 
tions spanning a range of no more than m adjacent lattice sites are retained; and (3) the 
localised LSUBm scheme, in which all multi-spin correlations over all distinct locales on the 
lattice defined by m or fewer contiguous sites are retained. Equation (j4]) shows that the 
ground-state energy at the stationary point has the simple form 



In this article we focus on the application of high-order (CCM) [22-24] to the spin-half, 
one-dimensional J1-J2 model. The CCM calculations were carried out to high-order using 
a computational code written by us [25]. The first CCM analyses [8,9] concentrated the 
Neel model state with neighbouring spins antiparallel [8] and then for a spiral model state 
[8,9] in which there is a pitch angle between neighbouring spins. It was assumed that 
translational symmetry was preserved in both cases, and these calculations led to good 
results for the ground-state energies of this system. In Refs. [10,11], we used a "doubled" 
unit cell including two neighbouring sites for this spin-half system on the linear chain at 
points (0,0,0) and (1,0,0) and a single Bravais vector (2,0,0) T to take into account the 
symmetry breaking. There are thus two distinct types of two-spin nearest-neighbor ket- 
state correlation coefficients, namely, those connecting the sites inside the unit cell and 



SH/SSj = ($|C7e- 5 #e 5 |$) = 0, V / ^ ; 
SH/SSj = ($\Se- s [H,Cf]e s \$) =0, V / ^ . 



(4) 
(5) 



E g = E g ({Sj}) = ($|e-^e s |$) 



(6) 
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those connectin g di fferent unit cells. These coefficients are denoted as <Sf and S%, and it is 



straightforward 
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llj to prove that the ground state at J2/J1 = \ is obtained exactly by 



setting SS; = 1 and all other coefficients equal to zero. Interestingly, a similar exact ground 
state at J2/J1 = \ may also be formed for the n.n.n. Neel model state and again S$ = 1 
and all other correlation coefficients are zero. 

We use three model states, namely, a spiral model state, and n.n. and n.n.n. Neel model 
states, shown in Fig. [TJ For all of these model states, we rotate the spin coordinates of the 
'up' spins so that notationally they become 'down' spins in these locally defined axes. For 
example, the relevant Hamiltonian in these rotated coordinates for the model state (a) is 
then given by 

H = E AA co <^) s t^ + \^os(e h] ) + l)(sr s J + sfsj) 

+-(008(0;, j) - l)(s+ S + + S r s j) + -sm(9 i!j )((sr s +) s * - s*( a r + (7) 

where 9^ is the difference between pitch angles at sites i and j. Furthermore, we see for 
model state (a) that #2i,2i+i = a, #2i+i,2i+2 = P, and 9 i)i+2 = a + f3 where i = 0, 1, 2, ■ ■ •. 
(We note also that J^+i = J\ and J^i+2 = J2) We remark that both model states (b) may 
be obtained from model state (a) by setting a = (3 = tt and model state (c) by a = it and 
(3 = 0. Model state (a) allows us to investigate the stability of the collinear n.n. and n.n.n. 
model states, as well as the special case a = (3 previously considered in Refs. [8,9]. 



III. RESULTS 

For J2/JI < |, the n.n. Neel state is found to yield the best results. For J2/JI > |, 
there were found to be two stable solutions with respect to the "pitch" angles a and (3, 
namely: the collinear n.n.n. Neel model state (a = it and (3 = or vice versa), which 
breaks the symmetry of the lattice; and, a spiral state with a = f3, which therefore does not 
break the translational symmetry of the lattice. We may plot the angle a(=(3) for the spiral 
solution, and this is shown in Fig. We see that the spiral solution is "lost" at around 
J2/J1 ~ 0.7, as noted previously in Ref. [9]. There is a direct and finite step-change to 
a = it at this point. Furthermore, we see slight oscillations in a for J2/J1 ~ 1- It is unclear 
why this should occur, although perhaps this behaviour may indicate some rapid change in 
the ground state and/or that the model state is becoming increasingly inappropriate. 
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FIG. 2: Results for the pitch angle a(=ft) versus Jij J\ for the spiral model state (a). 



We may investigate these solutions further by fixing the value of ft and then varying a, 
see Fig. El For the n.n.n. Neel model state, we may set ft = for all values of J%j J\ and 
we see from this figure that the minimum energy solution is given by a = ti in this case 
(i.e., identically model state (c)). Furthermore, we may set the value of ft to be given by 
the value shown in Fig. [2] (at a given level of LSUBm approximation and for a given value 
of J2/J1) in order to investigate the "spiral" case. We note that the stable minimal solution 
is given by a = ft, which preserves the translational symmetry of the lattice. Note that no 
other intermediate solutions, e.g., a kind of double spiral in which ft = and a 7^ ir were 
seen to occur. For J%j J\ > |, there was either the translational symmetry-breaking case of 
model state (c) or the spiral state with a = ft, which preserves translational symmetry. 

We remark that the n.n. Neel model state (b) (i.e., a = ft = tt) generally better produces 
for J2/J1 < \ than the n.n.n. Neel model state (c), and vice versa in the region J2/J1 > §■ 
Results for the ground-state energies as a function of J2/J1 are shown in Fig. HI Of the 
collinear model states (b) and (c), results for the n.n. Neel model state (b) are shown for 
J2/J1 < \ and those of the n.n.n. Neel model and (c) for J2/J1 > \- Exactly at the MG 
point J2I ' J\ = \, we find that these two states have equal energy. Furthermore, this solution 
is the CCM symmetry-broken solution (more is said of this below), which forms an exact ket 
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FIG. 3: Results for the CCM ground-state energy versus pitch angle a as a function of J2/J1 at 
the LSUB6 level of approximation. Results on the right of the figure are for (3 = (so are for the 
n.n.n. Neel model state (c)) and those on the left are for f3 given in Fig. [2] are so are for the spiral 
case. 

eigenstate at this point. For J2/J1 > |, we find that the a = tt and /3 = (i.e., model state 
(c)) is lowest until J 2 ML equal to 1.02, 1.69, 1.29, and 1.76 at the LSUB4 to LSUB10 levels 
of approximation respectively. (The non-monotonic behaviour of this crossing is probably 
due to a "double" odd/even effect, where the "doubled" comes from the range of the J 2 
bonds.) Above this point, the energies for the spiral model state (a) with a = (3 appear to 
be lower. Symmetry breaking might therefore persist for some region for J2/J1 > \- The 
differences between the energies of these two states above this crossover point are very small 
however, e.g., it is of order 10~ 4 at the LSUB10 approximation. The difference only becomes 
large for J2/J1 < J2I J\\c 2 ^ where the symmetry-breaking n.n.n. Neel collinear model state 
(c) is clearly much lower than the spiral solution for a — ft. Furthermore, it is possible 
that the crossover goes to higher and higher values of J2 with increasing level of LSUBm 
approximation level. However, we cannot rule out the possibility that another state of more 
"exotic" order that is not accessible from the model states used here might occur just above 
J2/J1 — \- If so then this state might exhibit lower ground-state energies. 

Henceforth, we concentrate on results for the collinear n.n. Neel model state (b) for 
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FIG. 4: Results for the CCM ground-state energy as a function of Jij J\ for the collinear n.n. Neel 
model state (b) (for J2/J1 < \) and n.n.n. Neel model state (c) (for J2/J1 > \) and the spiral 
model state (a) with a = (3. 

J2/J1 < \ an d the collinear n.n.n. Neel model and (c) only for J 2 / ' J\ > \. Results for the 
ket-state correlation coefficients for the n.n. and n.n.n. Neel model states (b) and (c) are 
shown as a function of Jij J\ in Fig. [5j We remark again that both model states (b) and 
(c) provide an exact solution at J2/J1 = \, namely, where S?; = and all other ket-state 
correlation coefficients are zero. Starting from J2/J1 = \ we were able to track this exact 
solution (as seen previously also in Refs. [10,11]) within a certain LSUBm approximation for 
the n.n. Neel model state to lower values of J2/J1 < \- It was also found previously [10,11] 
that the solution (i.e. S2 = S\) having the full translational symmetry [8,9] was the only 
solution below a critical point J2/Ji\a(< \)- Above this point the CCM solution for the 
n.n. correlation coefficient bifurcates, as shown in Fig. [5j CCM results for the ground-state 
for the symmetry-breaking solution to the CCM equations for J2/J1 > J2I J\\ci compared 
better to exact diagonalisation results in this regime than those results for the "usual n.n. 
('Neel-type') solution" that do not break translation symmetries (see Refs. [10,11] for more 
details.) However, as was also noted before [10,11], the value for J 2 / Ji\ ci is much too large 
(e.g., J2I J\\c x — 0.436 at the for LSUB14 level of approximation). This value ought to be 
compared to the value of J2/J1 = 0.2411(1) the model exhibits a transition to a two-fold 
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FIG. 5: Results for the CCM ket-state correlation coefficients at the LSUB12 level of approximation 
for the one-dimensional spin-half J1-J2 model using the collinear nearest-neighbour (n.n.) Neel 
model state and the collinear next-nearest-neighbour (n.n.n.) Neel model state. 

degenerate gapped dimerized phase. 

As we increase the level of approximation the value for J2I J\\ Cl does seem to decrease, 
albeit slowly. In Ref. [11], it was proposed that results for the predicted energy gap at 
a given level of LSUBm approximation are overestimated, with good results given only 
after extrapolation. This might artificially stabilise the "usual CCM solution" tat preserves 
translational symmetry and that is known to be gapless Hence Jij J\ | Cl might be grossly 
overestimated. However, it is fair to say that this proposition is mostly speculation. We 
believe however that the detection of any such "symmetry-breaking" solutions via the CCM 
is strong point of the method. 

We note again that the exact solution at J 2 / J\ = \ may be obtained for the n.n.n. Neel 
model state (c) by setting <5>f = 1 and all other coefficients equal to zero. This solution may 
be tracked for J 2 / 1 J\ > \ and the results for the n.n. ket-state correlation coefficients are also 
shown in Fig. |5j We remark again that this is the minimum energy solution compared to the 
"spiral" solution up to some possible cross-over value J%j Ji\ C2 . Importantly, the ket-state 
correlation coefficient S% is non-zero (and S% is zero) for all values of J2 > \ for this model 
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FIG. 6: Results for the CCM bra-state correlation coefficients at the LSUB12 level of approximation 
for the one-dimensional spin-half J1-J2 model using the collinear nearest-neighbour (n.n.) Neel 
model state and the collinear next- nearest-neighbour (n.n.n.) Neel model state. 

state, which suggests that the symmetry-broken ground state might persist for some range 
of J 2 /Ji(> |). 

The nearest-neighbor bra-state correlation coefficient at the LSUBm level of approxima- 
tion at J2/J1 = \ has S 2 = 1/4 with m > 4. This is shown in Fig. [6] for the LSUB12 
level of approximation. We find that the nearest-neighbor correlation coefficient diverges as 
J2I ' J\ J2I ' J\ C \ an( f this is also shown in Fig. [61 Again, the usual ('Neel-type') solution 
(<S 2 a = S b 2 ) is obtained for J 2 /Ji < J 2 /Ji\ ci - A gain, a similar solution at J2/J1 = n ma y 
also be be found for n.n.n. Neel model state (c), again for which S 2 = 0.25 and all other 
coefficients equal to zero. In the limit J2/J1 — > 00, the bra-state correlation coefficient S$ 
goes to zero. Once more, qualitatively similar results are observed at other levels of LSUBm 
approximation for the bra-state correlation coefficients as a function of J 2 / ' J\. 

The dimer order parameter is given by d = (s 2i ■ (s 2i+ i — s 2i _i)). We find this quantity 
for the n.n.n. Neel model state and the results for J2 > \ are shown in Fig. [71 Results 
for the n.n. Neel model state show an unphysical divergence at the at the critical point 
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FIG. 7: Results for the dimer order parameter for the one-dimensional spin-half J1-J2 model using 
the collinear next-nearest-neighbour (n.n.n.) Neel model state. 

J2I J\\c^ which is due to a divergence in the NCCM bra-state correlation coefficients as 
shown above. (This is uninteresting and so is not shown here.) We remark that there is a 
maximum in d near to J2/J1 ~ 0.6. For example, at the LSUB12 level of approximation, 
there is a maximum value of the dimer order parameter d of d = 0.790 at J2/ J\ = 0.58. This 
compares well to the results of DMRG in Refs. [5,9] by visual comparison. In particular, 
it was noted in Ref. [9] that the maximum in d occurs at J2/J1 = 0.58, in agreement with 
these new CCM results. The results for the n.n.n. Neel model state appear to decay to 
zero in the limit J2/J1 00, as expected. However, they do not seem to agree (at given 
LSUBm level of approximation) with the results of DMRG [5] that predict a decay of d 
given by d — 2.283exp(— I.622J2/J1). Furthermore (and somewhat troublingly) , the CCM 
results also show unphysical oscillations in d as one varies J2I J\- The reason for this is 
unclear, although it might hint at strong changes occurring in the in the ground state at 
this point and/or that the n.n.n. Neel model state is possibly a poor choice. However, it is 
probably still fair to say that the CCM results for the dimer order parameter broadly follow 
the same pattern as in Ref. [5,9], although quantitative agreement for large J2/J1 is very 
poor. Higher orders of approximation and/or extrapolation might solve this problem. 
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FIG. 8: Results for the spin-spin correlation function via the CCM for the one-dimensional spin-half 
J\—Ji model using the collinear next-nearest-neighbour (n.n.n.) Neel model state at the LSUB12 
level of approximation. 

Spin-spin correlation functions may be found for the n.n.n. Neel model state (c). At 
LSUB12 these results are found to be converged to roughly a separation of r m 10, shown 
in Fig. |H1 and at LSUB14 they are are found to be converged to roughly a separation of 
r w 12. Note however that convergence is also somewhat affected by the value of J2/J1, 
with convergence with increasing level of LSUBm levels of approximation being quicker 
at lower values of J2I J\- The correlation length e may be found [5] by fitting the form 
ax -0,5 exp(— x/e) cos(#x + A) to the spin-spin correlation functions for the n.n.n. Neel model 
state (c) and these results plotted as a function of J2/J1 in Fig. [9] at the LSUB10 to 
LSUB14 levels of approximation. By visual inspection only, these results at the LSUB14 
level of approximation do not appear to compare too badly to those of DMRG presented in 
Ref. [5]. Indeed, the correlation length does seem to increase with increasing J2I J\- This is 
an interesting result for a collinear model state of the type given by the n.n.n. Neel model 
state (c), although some caution should be exercised in interpreting these results because 
we again include results for results only up to r = 8 at LSUB10, r = 10 at LSUB12, and 
r = 12 at LSUB12. We see that the discrepancy between the predicted correlation lengths 
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FIG. 9: Results for the correlation length via the CCM for the one-dimensional spin-half J1-J2 
model using the collinear next-nearest-neighbour (n.n.n.) Neel model state at the LSUB10 to 
LSUB14 levels of approximation. 

between LSUB10, LSUB12, and LSUB14 becomes larger as J2/J1 increases. However, all 
levels of approximation indicate an increase in the correlation length with increasing J2/J1, 
as expected. 

The static structure function may be defined by 

1 r=+oo 

C(q) = ~ E (*(0) " *(r))e ir q (8) 

r=— 00 

However, again our localised LSUBm approximation will function correctly only for la re- 
stricted range of r. Results for the static structure function at LSUBM are shown in Fig. |9j 
We see that the peak of the structure function moves from q = ir at J%j J\ = 0.5 to q = n/2 
:or larger values of J^j ' J\- This is in reasonable agreement with previous results of Ref. 
9| that shows a similar shift in the peak, which was interpreted as indicating spiral order. 
However, the present results indicate that such a change can occur for a model state with 
collinear n.n.n. ordering also. By visual comparison, the heights of the peaks shown in Fig. 
[10] are in reasonable qualitative agreement with those shown in Fig. 4 of Ref. [9]. We remark 
that characteristic and strong oscillations are seen at the LSUB14 level of approximation 
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FIG. 10: Results for the static structure function C(q) via the CCM for the one-dimensional 
spin-half J1-J2 model using the collinear next-nearest-neighbour (n.n.n.) Neel model state at the 
LSUB14 level of approximation. 



for J2/J1 ~ 1- This might be the model state is a poor choice in this regime, although it 
is more likely that it is because we need to truncate the range r for the LSUBm scheme. 
Indeed, similar oscillations are seen in Fig. 4 of Ref. [9] and these are most probably due to 
finite-size effects. 



IV. CONCLUSIONS 



In this article, CCM results were presented for the spin-half ID J\-Ji model. We saw that 
the n.n. Neel model state provided good results for J2/J1 < \- In particular, a translational 
symmetry-breaking solution was shown to provide good results for the symmetry-broken 
phase for this model that exhibits a gap (see also Refs. [10,11]). However, as remarked 
previously elsewhere [10,11] the value for the onset of this phase J2I J\\c x = 0.2411 was 
shown to be overestimated by this approach. The stability of the (non-collinear) spiral 
and (collinear) n.n.n. Neel model states were examined by using a model state that could 
interpolate between the two model states (and incidentally also the n.n. Neel model state). 
It was found for J2/J1 > \ that the only stable solutions were for the case of the "non- 
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translational-symmetry-breaking" spiral and the (symmetry-breaking) n.n.n. Neel model 
state. All other intermediate cases seemed not to be stable. Interestingly, the ground-state 
energies for these two states were extremely close for large range of J^/Ji- However, the 
spiral case appeared to be slightly lower in the limit of large values of Jij J\. Even more 
startlingly, however, the energies for the collinear n.n.n. Neel model state became lower 
than for the spiral case in the regime \ ^ J21 3\ ~ 2. This hints at a possible "order-from- 
disorder" phenomena in this model in this regime. However, these results should be treated 
with caution because the "true" ground state in this regime might still be quite different to 
any of the states studied here. 

Results for the dimer order parameter d using the n.n.n. Neel model state showed a max- 
imum at d at about J 2 / 1 J\ ~ 0.58 (of d ~ 0.79 at the LSUB12 level of approximation), which 
was in good agreement with DMRG results of Ref . [9] . The decay of d with increasing J 2 / J\ 
did not agree with DMRG results and also the CCM results also demonstrated "oscillations" 
in d as one varied J2I J\- The n.n.n. nature of the J2 bond means that the length scale is 
effectively doubled and so we probably need to go to much higher levels of LSUBm approx- 
imation via massive parallel processing. In order to achieve good quantitative agreement 
for large J2/J1, we might also need to extrapolate our LSUBm results in the limit m — > 00 
However, it is probably still fair to say that the CCM results for the dimer order parameter 
d showed here followed a broadly similar pattern to results of DMRG of Refs. [5,9], 

The presence of collinear n.n.n. Neel ordering some region for J2/J1 > \ seems to 
contradict earlier results of other approximate method that indicate incommensurate spiral 
correlations commonly (see, e.g., Ref. [9]), i.e., possibly suggesting some form of "spiral" 
ground state. Spiral ordering is suggested because the classical spiral ground state must also 
show a similar change position of the peak in the (classical) structure function. However, 
we note that our results the static structure function for the n.n.n. Neel model (e.g., at the 
LSUB14 level of approximation in Fig. [TUj) state also showed a distinct peak that moved 
from q = n at J2/J1 — \ to q = n/2 as J2/J1 =— > 00. This was also seen in DMRG studies 
[9] of the quantum system. Furthermore, CCM results here predict that the correlation 
length becomes larger as one increases J2/J1 in the regime J2/J1 > \, again as was also 
seen in DMRG studies [5]. Note however that we used a fairly restricted range in r for the 
spin-spin correlation functions in order to form our results for the structure function etc. 
(only up to a separation of r « 10, 12 at LSUB14 for the structure function). This fact does 
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not bolster our confidence in our results, although this approach ought to be adequate in the 
regime just above Jij J\ = |, where the correlation length is known to be very short-ranged. 

In conclusion, the results presented in this article present an intriguing (though ultimately 
inconclusive) possibility that collinear n.n.n. Neel ordering might occur for J2/J1 > \- 
The ground-state energies for the n.n.n. Neel model state are lowest and results for the 
dimer order parameter, correlation length and structure function using this model state 
are (arguably) predicted correctly. If this result is in fact true then it has a number of 
implications for the interpretation of similar evidence (e.g., from structure functions) in 
other systems where the existence incommensurate "spiral" ordering is inferred. We suggest 
that other approximate methods (e.g., spin- wave theory and cumulant series expansions) 
might be applied to the ID J\—Ji model using a n.n.n. Neel state as a "starting point". We 
recommend also that even higher-order CCM be carried out using massive parallel processing 
in order to confirm the results presented here. 
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